Stéphane Clément 
Fernand Didier 
Dominique Proudhon 
Robert Rolland 
Patrick Soubeyrand 


TRAVAUX PRATIQUES DE 
MATHÉMATIQUES 
(CLASSES DES LYCÉES) 
TOME I 


Stéphane Clément 
E-mail : stephane.clement@ac-aix-marseille.fr 


Fernand Didier 
E-mail : didier@irem.univ-mrs.fr 


Dominique Proudhon 
E-mail : dproudhon@wanadoo.fr 


Robert Rolland 
E-mail : rolland@iml.univ-mrs.fr 
Url : http://iml.univ-mrs.fr/"rolland 


Patrick Soubeyrand 
E-mail : soubeyrand.patrick@free.fr, 


IREM d’Aïx-Marseille, Luminy Case 901, F13288 Marseille CEDEX 9 


26 Mai 2007 


TRAVAUX PRATIQUES DE 
MATHÉMATIQUES (CLASSES DES 
LYCÉES) 

TOME I 


Stéphane Clément, Fernand Didier, 
Dominique Proudhon, Robert Rolland, 
Patrick Soubeyrand 


AVERTISSEMENT 


Le but de cette brochure est de décrire quelques thèmes, d’où il est 
possible d'extraire des travaux pratiques de mathématiques aptes à faire 
pratiquer aux élèves des activités intéressantes dans le domaine des 
mathématiques, des applications des mathématiques et des interactions 
avec d’autres disciplines. On espère que ces activités pourront éclairer 
et enrichir les notions qui interviennent dans les programmes des classes 
des lycées. 


Par rapport au cours traditionnel qui fait un vaste balayage assez su- 
perficiel d’un domaine relativement large, le TP fournit un temps où 
l’élève peut entrer en profondeur dans un problème plus ponctuel, en 
décortiquer plus finement les mécanismes. Il peut s’apercevoir que dès 
qu’on s'intéresse à la mise en pratique des objets et des méthodes, les 
choses ne sont pas aussi simples ni parfois d’ailleurs aussi difficiles que 
ce qu'il aurait pu sembler à première vue. De nouveaux problèmes se 
posent, qui, sans la mise en pratique n’auraient pas pris un grand sens. 


Les TP tels que nous les concevons comportent plusieurs parties qui 
peuvent éventuellement ne pas être toutes présentes simultanément sur 
un même sujet : 

e Une présentation du problème avec son ou ses contextes scienti- 
fiques. 

e Une partie théorique et algorithmique sur les objets mathématiques 
utilisés et éventuellement les objets des autres disciplines impliqués 
dans le contexte scientifique. 

e Une mise en place des outils nécessaires à l’étude du problème et 
une production de résultats. 
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e Une interprétation et une discussion sur les résultats obtenus. 
Cette partie peut amener à un aller retour avec la partie précédente. 
La mise en place des outils nécessaires à l’étude du problème peut utili- 
ser un système de calcul et de calcul formel et comporter éventuellement 
une toute petite part de programmation. Nous avons choisi pour cette 
partie un logiciel libre de calcul formel performant : xcas. Mais ceci n’est 
pas exclusif : l’outil doit s’adapter aux problèmes qu’on à en vue et pas 
le contraire. Le tout est de disposer d’un outillage suffisamment com- 
mode pour ne pas détourner l’attention du travail principal au profit de 
l’apprentissage du fonctionnement d’un logiciel. 


Le rôle du professeur est bien évidemment primordial dans la mise en 
place d’un choix de TP intéressant. Et c’est un rôle difficile, car il de- 
mande une hauteur de vue sur des questions qui ne sont pas toujours cen- 
trales dans la pratique habituelle d’un enseignant de mathématiques. Il 
demande un travail important pour calibrer les sujets à une juste taille. Il 
demande enfin un investissement sur le choix des outils logiciels, adaptés 
à chaque TP, qu’il faut fournir aux élèves afin que ceux-ci puissent mener 
à bien leur travail. 

Contenus de la brochure : 

Les notes qui suivent décrivent des thèmes de travaux pratiques, 
ainsi que des développements scientifiques, remarques, compléments, 
implémentations, tests informatiques, nécessaires à notre avis aux pro- 
fesseurs pour cerner le sujet et pouvoir en fonction de sa classe proposer 
à partir de chacun de ces thèmes et de cet outillage un ou plusieurs TP 
adaptés. Il nous semble indispensable de fournir les tenants et abou- 
tissants des problèmes posés sans masquer l’environnement scientifique 
dont ils sont tirés. 

Cette brochure constitue donc un environnement de travail pour le 
professeur qui désire mettre en place une activité de travaux pratiques. Il 
ne faut donc pas considérer les sujets développés ici, comme des exemples 
à donner en l’état aux élèves, mais comme une aide pour les professeurs 
à la conception de sujets et d'activités. 

De même, la partie programmation est à adapter par un partage 
réfléchi entre ce qu’on donne tout construit aux élèves et ce qu’on leur 
demande de produire. Dans de nombreux cas, ceux-ci disposeront de 
tout le logiciel nécessaire au lancement des essais et des simulations 
qu'ils décideront compte tenu du problème posé. 
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CHAPITRE 1 


LISTE DE SUJETS 


1.1. Introduction 


On donne ici brutalement une liste de sujets dont on fournit ensuite, à 
l’usage des professeurs, le contexte scientifique. Ces sujets sont redonnés 
dans chaque chapitre. Il ne s’agit donc ici que d’une table récapitulative. 


1.2. Liste 


1.2.1. Lignes trigonométriques liées au pentagone régulier. — 
Trouver une méthode exacte pour vérifier par le calcul exact que : 


(2) 


) 4 


(On utilisera du calcul formel donnant une formule exacte et pas des 
valeurs approchées. En particulier on n’utilisera pas les fonctions trigo- 
nométriques de la machine). 


1.2.2. Lignes trigonométriques. — Construire avec un tableur une 
table trigonométrique donnant les sinus et cosinus des angles compris 
entre 0 et r/2 de 3 degrés en 3 degrés sans utiliser les fonctions trigo- 
nométriques. On partira des angles remarquables et on appliquera les 
formules trigonométriques classiques. 


1.2.3. Calcul approché du sinus de 1 degré. — Calculer le sinus de 
1 en utilisant les trois méthodes indiquées dans le chapitre correspondant. 


2 CHAPITRE 1. LISTE DE SUJETS 


1.2.4. Exploration d’une suite. — Explorer la suite dont le premier 
terme est uo = 0 est qui est défini par : 


1 
Un+1 — 3 (aus na a), 


où a = sin(r/60). 


1.2.5. Exploration d’une suite (1). — Il s’agit d'étudier pour la 
suite définie par : 


Un+1 — Fu) 
où 
F(x) = x(2 — ax), 


ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 
de convergence, l’influence de la valeur initiale. 


1.2.6. Exploration d’une suite (2). — Il s’agit d'étudier pour la 


suite définie par : 


Un+1 — Fu), 


F)= 3 (r+2), 


ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 


où 


de convergence, l’influence de la valeur initiale. 


1.2.7. Exploration d’une suite (3). — Il s’agit d'étudier pour la 
suite définie par : 


Un+1 — Fu), 
où 


2 
sr à 
FE (x) De ) 
2xæ —1 
ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 
de convergence, l’influence de la valeur initiale. 
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1.2.8. Exploration d’une suite (4). — On considère la suite de 
terme général : 
1 il 1 


Comparer expérimentalement S, et 7?/6. 
Faire de même avec la suite de terme général : 


1 
Sn RUE 
n 
1.2.9. Accélération de convergence. — On considère la suite de 
terme général : 
1 1 1 


Comparer expérimentalement S, et 7?/6. 


Essayer de déterminer expérimentalement a de telle sorte que la suite 
de terme général : 
2 TT A 
n 
converge le plus vite possible. 
Le nombre a étant maintenant fixé à la valeur prévue précédemment, 


essayer de déterminer expérimentalement b pour que la suite de terme 
général : 


a. b 
Sn — —, 
M M 
converge le plus vite possible. 
1.2.10. Développement décimal de période paire. — Le nombre 


1/7 a pour développement décimal 


1 
= = 0.142857. 
d 


Si on considère la première moitié de la partie périodique : 142, et qu’on 
l’ajoute à la deuxième moitié de la partie périodique 857, on obtient : 
999. Cela se produit-il dans d’autres cas (où la période est paire) ? 
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1.2.11. Nombres de Mersenne généralisés. — Un nombre de Mer- 
senne généralisé de base a est un nombre de la forme 

a 

a—1l” 


(a est un nombre > 1). Les nombres de Mersenne relèvent du cas où 
= 2; 


On va se placer dans le cas où a = 10. Un nombre de Mersenne de 
base 10 s’écrit sous la forme 11:--1 avec n chiffres 1. 


Peut on trouver des nombres de Mersenne à base 10 qui soient des 
nombres premiers ? 


1.2.12. Nombres de Mersenne généralisés (bis). — Quels sont les 
nombres qui divisent un nombre de Mersenne généralisé de base 10 (un 
nombre de la forme 111---1)? 


1.2.13. Nombres de Sophie Germain. — Peut-on trouver des 
nombres premiers q tels que p = 2q + 1 soit aussi un nombre premier 
(nombre de Sophie Germain) ? 

Dans ce cas quelle est la période du développement décimal de 1/p? 


1.2.14. Période maximale. — Existe-t-il des nombres premiers p tels 
que la période du développement décimal de 1/p soit égale à p—1? Dans 
un tel cas, soit 1 < k < p un nombre entier. Comparer le développement 
décimal de k/p et de 1/p. 


CHAPITRE 2 


CONSTRUCTION D’UNE TABLE 
TRIGONOMÉTRIQUE 


2.1. Présentation du problème 


2.1.1. Les grandes lignes. — Le but de ce travail est de tabuler 
sin(x) et cos(x) de degré en degré entre 0° et 360°, c’est-à-dire pour les 


angles de la forme : 
kT 


TR — Sr: 
180 
où k est entier et vérifie 0 < k < 360 (k est la mesure de l’angle en degré). 
Des considérations géométriques nous permettent de connaître les sinus 
et cosinus d’un certain nombre d’angles particuliers 0°,18°, 30°,45° : 


FIGURE 1. Lignes trigonométriques classiques 


Compte tenu des formules classiques : 


sin(x + 2kr) = sin(x), 
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T 
cos(s — x) = sin(x), 
on se restreint aux angles x tels que : 
T 
O<x<—, 
RE: 

c’est à dire : 

0 < k < 45. 
2.1.2. Disgression géométrique. — Les valeurs de la table 1 


découlent directement de considérations géométriques classiques : pour 
l’angle 0 il n’y à évidemment pas de problèmes. Le sinus et cosinus de 
m/6 se calculent en considérant un triangle équilatéral, pour 7/4 on 
considère un triangle rectangle isocèle. Le cas de 7/10 est un peu plus 
compliqué et utilise le pentagone régulier ou si on préfère les racines 
cinquièmes de l’unité. Notons donc : 


2im 
DES, 
Les racines cinquième de l’unité sont alors : 


4 
1, a, a, a*, a | 


Calculons la somme de ces racines : 


S=1+a+ a? ++. 


Cette somme en vertu de la formule de sommation des premiers termes 
d’une progresson géométrique vaut : 
5 


ne 


1—-a° 


et comme a° = 1 on obtient : 


Remarquons que si on note : 
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alors les racines cinquième de l’unité s’écrivent aussi : 


1,, 8, B, a, 


et donc : 


S=1+(a+a)+(8+8)=0, 


ou encore en utilisant les formules d’Euler : 


27 Ar 
1+2 — 2 —.} = 40; 
+ co (2) + co (TE) 
4 2 
COS (=) — 2 cos? (=) — 1, 
et en conséquence : 


2 2 
4 cos? (=) + 2 cos (2) — 1 =0. 


Ainsi cos(27/5) est la solution positive de l’équation : 


4x? + 92% — 1 =0. 


Mais : 


On trouve : 


On a alors aussi : 


c’est-à-dire : 


2.2. La table trigonométrique de 3° en 3° 


À partir des valeurs fournies par la table 1 et des formules trigo- 
nométrique : 


(1) sin(a — b) — sin(a) cos(b) — sin(b) cos(a), 
(2) cos(a — b) = cos(a) cos(b) + sin(a) sin(b), 


on peut calculer successivement les lignes trigonométriques des angles 
15° (45° — 30°) et 12° (30° — 18°) puis 3° (18° — 15°). 
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Ainsi maintenant grâce aux formules : 


(3) sin(a + b) = sin(a) cos(b) + sin(b) cos(a), 


(4) cos(a + b) = cos(a) cos(b) — sin(a) sin(b), 


appliquées avec b — 3r/180, on peut calculer formellement la table des 
sinus et des cosinus de 3° en 3°. On peut obtenir des formules “exactes” 
écrites avec des radicaux. Ces formules sont assez lourdes, mais peuvent 
être trouvées effectivement avec un logiciel de calcul formel. 


On peut bien entendu traduire ces formules exactes avec radicaux en 
des formules approchées par des écritures décimales. On peut d’ailleurs 
mener les calculs dès le début avec des approximations décimales. 


2.3. La table trigonométrique de degré en degré 


Puisqu’on dispose de la table de 3° en 3° il suffit de calculer le sin et 
le cos pour un angle d’un degré et d’utiliser les formules d’addition et 
de soustraction pour avoir tout ce qu’on cherche. Notons u le sinus de 
l’angle r/180 et v son cosinus. On connaît déjà les sinus et cosinus des 
angles de la forme 3kr/180 où 0 < k < 15. On calcule alors pour k = 0 
jusqu’à k = 14: 


(5) sin En oS = sin — U + u cos si 
180 : 180 180 / ? 


(6) cos LE 608 ul v — usin cet 
180 _ 180 180 / ” 
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2.3.1. Calcul du sinus et du cosinus de 1° par réinjection. — 
On part de la formule : 


(9) sin(3x) = 3sin(x) — 4sin*(x) 
qu’on va appliquer ici avec x = 7/180. On cherche alors à résoudre cette 
équation connaissant la valeur a = sin(3x). L’équation se ramène à : 
1 

sin(x) = 3 sin?(x) + a). 

Donc si on pose : 
1 
J(u) = (au + a), 
le nombre sin(r/180) est solution de l’équation : 
fu) = u, 

c’est-à-dire est un point fixe de f. Pour calculer cette solution on part 
d’une valeur approchée de la solution wo, par exemple wo = 0, et on calcule 
ui — f(uo). La valeur u1 est réinjectée dans le second membre pour 


calculer une nouvelle approximation u2 — f(u1), puis plus généralement, 
par récurrence : 


Ua = Jui} 


Théorème 2.3.1. — La suite de terme général u, converge vers la va- 
leur cherchée sin(r/180). 


Démonstration. — La fonction f a pour dérivée : 
f'{(u) = du”, 


donc elle est croissante. 
On à par ailleurs : 


et aussi : 
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donc : : : 

10 
On conclut que l’image par f de l’intervalle [0, a/2] est incluse dans l’in- 
tervalle [0, a/2], autrement dit, restreinte à l'intervalle [0, a/2], f est une 
application de [0, a/2] dans lui-même et cette application est croissante. 
Comme uo = 0, u1 = f(uo) > Uo, on à aussi f(u1) > f(uo), c’est à dire 
que u2 > u et par récurrence un > Un_1. La suite de terme général uy 
est croissante, elle est majorée par a/2, elle converge donc vers une limite 
b. Cette limite vérifie b = f(b). L'étude de la fonction g(u) = u — f(u) 
sur l’intervalle [0, a/2] nous montre que l’équation g(u) = 0 à au plus une 
solution dans cet intervalle (en effet g(u) est strictement croissante sur 
cet intervalle). Donc le b trouvé comme limite de la suite de terme général 
un est l’unique solution de l’équation u = f(u) sur l'intervalle [0, a/2]. 
Comme on sait que sin(r/180), qui est dans l’intervalle en question, est 
solution de cette équation, pas de doute, b = sin(r/180). Ü 


On peut aussi évaluer la qualité de l’approximation. Sur l'intervalle 
Drap 
LP'(x)| < a°. 
On à alors en utilisant l'inégalité des accroissements finis et le le fait que 
D f{ble 
b— un] = 1f(0) — f(un-1)| € ab — uni], 


et en réitérant le procédé, 
b— un] < (a°)*1b — wo], 


et en majorant |[b — wo| par la longueur de l'intervalle : 


2n+1 
b—u,| < 


2.3.2. Calcul du sinus et du cosinus de 1° par une approximation 
polynomiale. — On va prendre ici pour valeur approchée du sinus au 
voisinage de 0 le polynôme : 


P(x) in 


Compte tenu du problème posé on prendra x > 0. 


2.3. LA TABLE TRIGONOMÉTRIQUE DE DEGRÉ EN DEGRÉ 11 


Pour évaluer la qualité de cette approximation nous allons partir de 
l'inégalité triviale 
cos(x) < 1. 


Par conservation des inégalités par intégration sur un intervalle {a, b] où 


b > a, on a donc : 
l cos(t) dt < | dé, 
0 0 


sin(t) <r: 


c’est-à-dire : 


En réappliquant cette méthode on obtient maintenant : 


ou encore : 


et en conséquence l’encadrement suivant : 


_ 
1 — Se < cos(x) < 1. 


On obtient alors successivement, toujours par intégration les encadre- 
ments suivants : 
Re 
Fee <sin(r) 


ia L# LE 
le See 
< cos(x) < 5 +52 
3 He - GE 
Ne pr 
æ ES en & + 125 
Ceci montre en particulier que : 
15 
0 < si — P < —. 
< sin(x) Gohe 5 


Appliqué au cas où x — r/180, on obtient : 


T T 
UE -P() <1 10711. 
0<sn( =) 180) < 1-35 10 
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2.4. Étude expérimentale par informatique 


Commençons par une petite vérification : on vérifie par le calcul formel 
que : 
V5 —1 
4 
est bien égal à sin(r/10) c’est-à-dire à cos(27/5). 


On entre interactivement les 4 lignes de commandes suivantes au 
prompt de xcas : 


a:=(sqrt(5)-1)/4 
b:=sqrt(1-a°2) 
z:=a+(I)xb 
simplify(z°5) 

La réponse est 1. Ceci vérifie le résultat. En effet, a et b sont donc 
respectivement le sinus et le cosinus d’un angle 0 tel que le nombre com- 
plexe de module 1 et d’argument © soit une racine cinquième de l’unité. 
Comme a et b sont tous les deux positifs strictement, l’argument 0 vaut 
27/5. 


On construit maintenant par programme la table trigonométrique. On 
va travailler dans deux listes ayant de 46 éléments. La première liste est 
constituée des 46 valeurs du sinus pour les angles de 0° à 45°, alors que 
la deuxième liste contient les cosinus de ces mêmes angles. Attention, 
en émulation maple, les listes partent de l’indice 1. Ceci introduira un 
décalage. Ainsi par exemple, sin(r/6) se trouvera dans la case numérotée 
31 de la première liste. 


Dans un premier temps on travaille en calcul exact, ce qui n’introduit 
que des racines carrées. De ce fait on se contente du calcul des lignes trigo- 


nométriques des angles de pris de 3° en 3°. Seules les cases 1,4,7,---,3k+ 
1,---,46 de la première et de la deuxième ligne contiennent des données 
intéressantes. 


## Calcul exact 
finitialise:=proc() 
local u,v,s,c,i; 

u:=0; 
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for i from 1 to 46 


u:=u,0; 
v:=v,0; 
od; 

s:= [ul]; 
c:=[v]; 
s [1] :=0; 
ELLES 
s[191:=(sqrt(5)-1)/4; 
c[19]:=(sqrt(10+2*sqrt(5)))/4; 
s[31]:=1/2; 
c[31]:=(sqrt(3))/2; 
s [46] :=(sqrt(2))/2; 


c[46]:=(sqrt(2))/2; 

s[13]:= (s[31]*xc[19]-s[19]*xc[31]); 
c[13]:= (c[31]*xc[19]+s[31]*xs[19]) ; 
s[16]1:= (s[46]xc[31]-s[31]*xc[46]) ; 
c[161:= (c[461xc[31]+s[46]*xs[31]) ; 


s[41:= (s[16]1*xc[13]-s[13]1xc[16]) ; 
cf4l:= (c[16]1*xc[13]+s[16]1xs[13]) ; 
for i from 1 to 14 


do 
if s[3*i+1]=0 
then 
s[3*xi+1]:= (s[3*i-2]xc[4]+s[4]xc[3*xi-2]) ; 
cl3xi+i]:= (c[3*i-2]xc[4]-s[4]xs[3*xi-2]) ; 
fi; 
od; 
RETURN([s,c]); 


end; 


Maintenant, on donne les valeurs en nombres réels. Le calcul du sinus 
de 1° est fait par itération ainsi que décrit dans le paragraphe 2.3.1 


## Calcul approché 


14 
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initialise:=proc() 


local w,t,s,c,sn,cn,i; 


w:=0.0; 
t:=0.0; 
for i from 1 to 46 

do 

wWw:=wW,0.0; 
t:=t,0.0; 

od; 
sn:=[w] ; 
cn:=[t]l; 
sn[1]:=0.0; 
cn[1]l:=1.0; 
sn([19]:=evalf((sqrt(5)-1)/4); 
cn[19]:=evalf((sqrt(10+2xsqrt(5)))/4); 
sn[31]:=0.5; 
cn[31] :=evalf((sqrt(3))/2); 
sn [46] :=evalf ((sqrt(2))/2); 
cn [46] :=evalf ((sqrt(2))/2); 
sn[13]:= (sn[31]x*xcn[19]-sn[19]*cn[31]); 
cnl13]:= (cnl31l*xcn[19]+sn[31]x*xsn[19]) ; 
sn[16]:= (sn[46]*xcn[31]-sn[31]*xcn[46l) ; 
cn[16]:= (cn[46]*xcn[31]+sn[46]x*sn[31]); 
sn[4]:= (sn[16]*xcn[13]-sn[13]*xcn[16]l) ; 
cnl4]:= (cnl16]*xcn[13]+sn[16]x*xsn[13]l) ; 
for i from 1 to 14 


do 
if sn[3*i+1]=0 
then 
sn[3*xi+1]:= (sn[3*xi-2]*xcn[4]+sn[4]l*cn[3*xi-2]); 
cn[3*i+1]:= (cn[3*xi-2]*+xcn[4]-sn[4]x*xsn[3*xi-2]) ; 
fi; 


od; 


2.4. ÉTUDE EXPÉRIMENTALE PAR INFORMATIQUE 


RETURN([sn, cn] ); 
end; 


tables:=proc(iter) 
local i,u,v,a,t,sn,cn; 
t:=initialise(); 


a:=t [1] [4]; 

u:=0.0; 

for i from 1 to iter 
do 

u:=(4x(u°3)+a)/3; 

od; 

v:=sqrt(1-u"2) ; 

sn:=t{1]l; 

cn:=t{[2]l; 

for i from O to 14 
do 


sn[3*i+2]:=sn[3*i+1]*v+tu*xcn[3*i+1]; 
cn[3*i+2]:=cn[3*i+1]*xv-u*xsn[3*i+1]; 
sn[3*x(i+1)]:=sn[3*x(i+1)+1]xv-uxcn[3*(i+1)+1]; 
cn[3*x(i+1)]:=cn[3*x(i+1)+1]xv+uxsn[3*(i+1)+1]; 
od; 
RETURN([sn,cn]) ; 
end; 


## Exécutions 


On calcule par l’approximation polynômiale 
le sinus de un degré: 


v:=evalf(pi/180) 
0.0174532925199 


w:=v-(v" 3)/6 
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0.0174524064238 


On regarde ce qu’il y a dans la table 
pour le sinus de un degré: 


t:=tables(100) 
t[1] [2] 
0.0174524064373 


Si on rajoute à la valeur donnée par l’approximation 
polynômiale de degré 3 le terme suivant (en w puissance 5) on trouve 
ce qu’il y a dans la table. 


w+1.35*107(-11) 
0.0174524064373 


Jusqu’à quel point on peut 
diminuer le nombre d’itérations : 


pour 4 itération c’est encore bon! 
t:=tables(4) 

€ [1] [2] 

0.0174524064373 


pour trois itérations on perd de la qualité 
t:=tables(3) 

+ [1] [2] 

0.0174524064268 


Remarque : on peut aussi résoudre l’équation 
3u — Au — 0.52335956242943835e — 1 


avec xCas : 


solve(3*u-4*u"3-0.052335956242943838 , u) 
Réponse : 0.017452406437283513 
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2.5. Quelques sujets sur ce thème 
2.5.1. Lignes trigonométriques liées au pentagone régulier. — 
Vérifier que : 
27 V5-—1 
COS | — | = ———. 
6! 4 


2.5.2. Lignes trigonométriques. — Construire avec un tableur une 
table trigonométrique donnant les sinus et cosinus des angles compris 
entre 0 et r/2 de 3 degrés en 3 degrés sans utiliser les fonctions trigo- 
nométriques. On partira des angles remarquables et on appliquera les 
formules trigonométriques classiques. 


2.5.3. Calcul approché du sinus de 1 degré. — Calculer le sinus 
de 1 en utilisant les trois méthodes indiquées. 


2.5.4. Exploration d’une suite. — Explorer la suite dont le premier 
terme est uo = 0 est qui est défini par : 


1 
UnH1 — 3 (dur + a), 


où a = sin(r/60). 


CHAPITRE 3 


CALCUL DE L’'INVERSE D'UN 
NOMBRE 


3.1. Présentation du problème 


Soit a un nombre réel que nous supposerons > 0 dont nous vou- 
lons calculer l'inverse 1/a. La méthode que nous présentons, et qui est 
implémentée dans certaines calculatrices et certains langages de program- 
mation, est très performante. C’est la méthode de Newton, développée 
pour le cas particulier qui nous préoccupe. 


3.2. La méthode 
Introduisons la fonction : 
F(x) = x(2 — ax). 
Cherchons les solutions de l’équation : 
(10) Gr (s) 
c’est-à-dire de : 
(11) x(1 — ax) =0 


dont les solutions sont 0 et 1/a. L'idée est de partir d’une valeur ap- 
prochée uo > 0 de 1/a de calculer u; = F(uo), de réinjecter u; dans la par- 
tie droite de l’équation 10 pour obtenir u2 = F(u;) et plus généralement 
par récurrence u, = F(u»_1). 

Notons J l’intervalle ouvert ]0,2/al. 
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Théorème 3.2.1. — Soit uo € I. Alors la suite définie par le premier 
terme wo et la relation de récurrence pour n > 1 
Un — F(un_1), 


converge vers 1/a. 


Démonstration. — La dérivée de F(x) sur l’intervalle Z est donnée par : 
F'(x) = 2(1 — ax). 


Sur l'intervalle 7, la fonction est positive et atteint son maximum en x = 
1/a. En ce point la valeur de la fonction est 1/a. On peut donc conclure 
que l’image de l'intervalle Z est incluse dans 7. Remarquons aussi que sur 
l'intervalle ]0, 1/a] la fonction F(x) est strictement croissante. Si wo € I 
alors 0 < F(uo) < 1/a. Posons ensuite u1 = F(uo) et par récurrence 
Un — F(un-1). Tous les termes u1,u2,--: sont donc dans l'intervalle 
10, 1/a]. De ce fait, au, < 1, si bien que (2 — au,) > 1 et donc 


tn = CG) > 


La suite (u,),>1 est croissante, majorée par 1/a, donc converge vers une 
limite 0 < b < 1/a qui vérifie b — b(2 — ab). Cette limite est donc le point 
1/a de F. Hal 


On peut étudier aussi la qualité de l’approximation obtenue. Pour cela 
on suppose qu’on parte d’une valeur w, telle que 


1 
0 <'uo < —. 
(a 


Dans ces conditions on sait que toutes les valeurs u, sont aussi dans cet 
intervalle. Alors on a successivement : 


1— au, = 1— auy_1(2 — aus_1), 
l— au, = 1 — 2aus_1 + Que ni 


lou, = (1-01), 


On obtient donc : 
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Posons à = 1 — aus. Compte tenu des conditions sur la valeur initiale wo 
le nombre à vérifie 0 < à < 1. La suite u, qui vérifie : 


a donc une convergence quadratique vers 1/a. 


3.3. Étude expérimentale par informatique 


3.3.1. Le choix d’un point de départ. — On commence par écrire 
a sous la forme : 

a = C 10, 
où 1 < c < 10 et où k est un entier relatif. Ceci permet de se ramener à 
calculer l’inverse d’un nombre compris entre 1 et 10. Ceci étant fait, on 
peut supposer désormais que 


1<a<10. 


Choisissons u, de la façon suivante : 

(1) si 1 <a < 2 alors üo = 0.5, 

(2) si 2 < a < 3 alors üo = 0.3, 

(SL Sn Dralore:us = 0; 

(4) si 5 < a < 10 alors u = 0.1. 
Ce choix respecte la condition u9 < 1/a, et de plus aus > 0.5. Donc 
1 — aug < 0.5. Ces choix étant faits on obtient puisque 1/a < 1 et 


a <1/2: 
IN 
< | — ; 


3.3.2. Programmation et calculs. — 


Ë 
a 


normalisation:=proc(x) 


local k,a; 

k:=0; 

a:=X; 

if (a>=1) 
then 


while (a>=10) 
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a:=a/10; 
k:=k+1; 


fi; 
RETURN([a,k]) ; 
end; 


F:=proc(x,a) 
RETURN (xx (2-axx) ) ; 
end; 


unsur:=proc(x,iter) 
local v,a,i,y,u; 
v:=normalisation(x) ; 
a:=v[1]; 
u:=0.5; 
if (a>=2) then u:=0.3; fi; 
if (a>=3) then u:=0.2; fi; 
if (a>=5) then u:=0.1; fi; 
for i from 1 to iter 
do 
u:=F(u, a) ; 
od; 
RETURN (ux(10°7(-v[21))); 
end; 
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3.4. Pourquoi la convergence est si rapide 


On peut déjà compte tenu de la situation penser que la convergence va 
être rapide. En effet compte tenu de l’inégalité des acroissements finis : 


Un — : — f(un_1) sf (:) < RU = U 
a a a 


où k est un majorant de F’(x) sur l’intervalle considéré. Donc : 


1 


+ 4 UQ — — 
a 


de 
a 


nl 


si bien que si k < 1 la suite converge. Mais ici, la dérivée est nulle au 
point fixe, donc on peut s’attendre à une amélioration de la rapidité 
d’approximation lorsqu'on se rapproche du point fixe. On à vu en fait 
que l’approximation est quadratique. 


O u0 ul 1 


FIGURE 1. Approximation de 1/a 
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3.5. Généralisons 


Soit f(x) une fonction sur un intervalle Z = [a,b] et ayant sur cet 

intervalle un seul zéro c. On suppose que f'(x) ne s’annule pas sur J. 

F()=5- F0) 

f(x) 
Alors c est le seul point fixe de F(x) sur l'intervalle J. 
1! 
sa =toie) 
f(x} 


Donc f'(c) = 0. On se retrouve dans la situation favorable précédente. 


A(x) 
F(z) — F(c) = F(x) —c Fa) 
Comme F(c) — c = 0 on conclut que A(c) = 0. Comme F"(c) = 0 on 
conclut que A’(c) = 0. Si par exemple A(x) est un polynôme ça veut dire 
qu'il s'écrit : 
A(x) = (x — a) B(x), 

où B(x) est un polynôme. 

Exemple : 

Prenons : 

f(n = -r-1 

et regardons ce qu'il se passe sur l'intervalle [1,2]. On est alors amené à 
introduire : 


ge LT 
ire | 
GE 


3.6. Quelques sujets sur ce thème 
3.6.1. Exploration d’une suite (1). — Il s’agit d'étudier pour la 
suite définie par : 
Un+1 = F(un), 
où 
F(x) = x(2 — ax), 
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ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 
de convergence, l’influence de la valeur initiale. 


3.6.2. Exploration d’une suite (2). — Il s’agit d'étudier pour la 
suite définie par : 
Un+1 — Fu); 


F)= 3 (r+2), 


ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 


où 


de convergence, l’influence de la valeur initiale. 


3.6.3. Exploration d’une suite (3). — Il s’agit d'étudier pour la 
suite définie par : 
Un+1 — Fu 
où 
x? +1 
2x — 1 
ainsi que par la donnée d’une valeur initiale wo, la convergence, la vitesse 


Flx)= 


de convergence, l’influence de la valeur initiale. 


CHAPITRE 4 


STATISTIQUES SUR LES PGCD 


4.1. Présentation du problème 


On tire au sort deux nombres entiers naturels a et b ayant au plus 
1024 bits, c’est-à-dire des nombres < 21%, De tels nombres ont approxi- 
mativement 300 chiffres décimaux. On suppose que chaque entier n tel 
que : 

O<n<20%42] 


a la même chance que les autres d’être tiré. On se fixe un entier h. Quelle 
est la probabilité P(h) pour que : 


pgcd(a, b) > h. 


4.2. Expérimentation par informatique 


toutalea:=proc(x) 
local n,s; 
s:=0; 
for n from 1 to x 
do 
s:=rand(2)+2x*s; 
od; 
RETURN (s) ; 
end; 


aleagcd:=proc(x) 
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local a,b,d; 
a:=toutalea(x) ; 
b:=toutalea(x) ; 


d:=gcd(a,b); 
RETURN (d) ; 
end; 


## x est le nombre de bits maximum choisi 

## niv est le nombre choisi si pgcd(a,b) >=niv on le note 

## n nombre d’essais 

## en sortie la proba (approchée) pour que le gcd de 2 nombres 
## de x bits tous tirés au sort soit >= niv 


gcdstat:=proc(x,niv,n,alea) 
local i,total,a,b,d; 
srand(alea) ; 
total:=0; 
for i from i ton 
do 
a:=toutalea(x) ; 
b:=toutalea(x) ; 
d:=gcd(a,b) ; 
if d>= niv 
then 
total:=total+1; 
fi; 
od; 
RETURN (evalf (total/n)) ; 
end; 


gcdstati:=proc(x,niv,n) 
local i,total,a,b,d; 


4.2. EXPÉRIMENTATION PAR INFORMATIQUE 


total:=0; 
for i from 1 to n 


a:=toutalea(x) ; 
b:=toutalea(x) ; 
d:=gcd(a,b) ; 
if d>= niv 
then 
total:=total+1; 
ÉT: 
od; 
RETURN (evalf (total/n)) ; 
end; 


lance:=proc() 

local mi,m2,m3,m4,m5,m6,m/,m8,m9; 

srand(427) ; 

= gcdstat1(1024,2,10000) ; 

m2:= gcdstat1(1024,3,10000) ; 

m3:= gcdstat1(1024,4,10000); 

md:= gcdstat1(1024,5,10000) ; 

m5:= gcdstat1(1024,6,10000) ; 

m6:= gcdstat1(1024,7,10000); 

m7:= gcdstat1(1024,8,10000) ; 

m8:= gcdstat1(1024,9,10000) ; 

m9:= gcdstat1(1024,10, 10000) ; 

RETURN([mi,m2,m3,m4,m5,m6,m7 ,m8 ,m9] ) ; 
end; 


Une valeur approchée théorique pour P(h) est : 


29 


Il semble que la probabilité P(h) soit peu sensible au nombre de bits 


fixé pour les tailles maximales de a et b, dans la mesure où on reste dans 


des grandes tailles. 
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4.8. Un peu de théorie 
4.3.1. Les notations. — Soit N un entier naturel. On pose : 
TEL SEAUNE 
On définit alors les ensembles suivants : 
E={(ab)el?|anb=1}, 
A={(ab)el*|b<aetaAb=1}, 
B={(a,b)e |a<betaAb=1}. 


On remarque que le couple (1,1) est le seul couple de 7?, de nombres 
premiers entre eux qui ne soit ni dans À ni dans B. Avec les notations 
précédentes, on dispose des relations suivantes : 


ANB =, 
E = AUBU#{(1,1}, 
HA=%P, 
et donc : 
HE =2#A +1. 


Il convient de calculer le nombre d’éléments de À. Pour cela on introduit 
pour tout a fixé : 


A; ={bEel|b<aetaAb=1}. 
Alors : 
A={| JA x {a} 


a€cl 
De plus les ensembles 4, x {a} sont disjoints. Pour tout a > 2, on a : 


#Aa X {a} 27 # Aa + o(a), 
où ® est la fonction indicatrice d’Euler, et clairement À; = Ÿ. Donc : 
#A= Ù é(a)-1. 
1<a<N 


On conclut que : 


#E=2 Ÿ g(a)-1. 


1<a<N 
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4.8.2. Probabilité de tirer un couple de nombres premiers entre 
eux. — Mettons sur l'intervalle 7 la probabilité équirépartie. On effec- 
tue l’expérience suivante : on tire au sort un nombre a € I et on tire au 
sort un nombre b € 1, ces deux tirages étant faits indépendamment. On 
cherche la probabilité p pour que a et b soient premiers entre eux. On a 
bien entendu : 


HE 
P — "N2° 
ce qui donne : 
_2 D 1<a<N (a) —1 
RS CO 
On sait que 
D icon g(a) 
N2 


est à peu près égal à 3/7°. On pourra voir que par exemple pour N > 10 
la probabilité p peut être prise raisonnablement pour : 
au 
Pp— 2 
4.3.3. Probabilité q pour que pgcd(a,b) > h. — Prenons N de 
l'ordre de 2° par exemple. 
Posons : 
Gy = {(a,b) € L° | pgcd(a, b) = k}. 
Il est clair que : 
Gr = {(a,b) € L° | a = k1k,b = kok, pgcd(k, ka) = 1}. 
Notons aussi : 


He = {(kiske) € F° | CRD ki A ko = 1}. 


Z 


Alors d’après les définitions on voit que : 


#Gr = #Hx. 
Mais on à approximativement : 
6 N? 


Si on introduit maintenant : 


G = {(a,b) € I? | pgcd(a,b) > h}, 
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alors : 5 | 
— Eu ERP = 2 — 
46= Lime dm y À 
R£k<N R<K<N 
Et la probabilité q vaut approximativement : 
6 1 
qT— 72 7.2? 
R<k<n 


ce qui est encore approximativement égal à : 
6 1 
mr? hk 


(évaluation asymptotique du reste de la série de terme général 1/k?). En 


conséquence on peut s’attendre à une probabilité approchée : 


4.4. Quelques sujets sur ce thème 


4.4.1. Concentration en nombres premiers. — Par une évaluation 


statistique essayer de prévoir la concentration en nombre premiers dans 


les entiers. 


4.4.2. Probabilité pour que deux nombres soient premiers entre 


eux. — Par des tests statistiques, conjecturer un résultat sur la pro- 


babilité pour que deux nombres entiers tirés au hasard soient premiers 


entre eux. 


4.4.3. Le pgcd de deux entiers. — Par des tests statistiques, conjec- 
turer un résultat sur la probabilité pour que deux nombres entiers pris 


au hasard aient un pgcd qui soit > d. 


CHAPITRE 5 


LIMITE D’UNE SUITE 


5.1. Présentation du problème 


Il s’agit d'étudier la convergence de la suite (5,),:, de terme général : 
1 1 1 
R=gtat te 
Puis, ayant établi la convergence de cette suite vers une limite S, nous 
essaierons d'évaluer la rapidité de la convergence en évaluant S — $,. 


5.2. Convergence de la suite 


La suite (5,),>, est croissante. En effet : 


1 
Snt1 — On = Gap 
donc : 
DR One 
Considérons la courbe (C) d’équation : 
#1 
Y — me. 


La somme S, s’interprête comme la somme des aires des triangles ha- 
churés sur la figure 1. En considérant par ailleurs l’aire sous la courbe, 
on en déduit que : 

Fdz 


Sn <1+ —. 
x 


1 
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En conséquence : 
1 
ne dise 
n 


La suite (S,),, étant croissante et majorée converge. Elle converge vers 


un nombre S < 2. On peut montrer que S — …. ce que nous admettrons 
ici. 


FIGURE 1. Majoration 


5.3. Évaluation de S — S, 


Posons : 
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L’entier n étant fixé, pour tout entier m tel que m > n +1, posons : 


1 L 


Re RES 
Ge Ne 


qu’on peut encore écrire : 
Res S or = Sr 


L’entier n étant fixé, la suite (Pie converge donc vers R, = S—S,. 


On obtient un encadrement de À, par une majoration du style de celle 
suggérée par la figure 1 et par une minoration telle que celle suggérée par 
la figure 2. 


1/4 


n n+]l 


FIGURE 2. Minoration 
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1 1 1 1 
Nerf Fa rs 
MEL, AE R CEE R 
et donc : 
res 
n+l on 
Cet encadrement nous permet d’obtenir l’encadrement suivant : 
Vi pele wrr 
TT on n+1. nr 


5.4. Amélioration de l'évaluation de S — 5, 


On va évaluer plus finement la différence entre l’aire sous la courbe et 
l’aire des rectangles construits sur la figure 1. Pour cela notons : 


et : 
td 
en) _ 
D 
Alors : 
is Re ÿ PES 
k=n+1 
et en passant à la limite sur m : 
/ R I S A 
——R, = lim 
nm m— 00 F 
k=n+1 
Or : 
1 1 1 1 


EI &k KR K(Kk-1l) 
ce qui nous permet de donner l’encadrement suivant : 
1 1 


UN Ve nes 
REED 
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En faisant cette fois-ci intervenir la courbe d’équation : 
1 
YU — 23° 
et en raisonnant sur des aires bien choisies comme dans le paragraphe 
précédent on obtient successivement : 


1/1 1 .—. jL 1 1 
2\R (+1) = °=2\(K-27 (k=-1ÿ/’ 
et par sommation et simplification : 
1 1 1 - d 1 1 
A PR D nn 2 à, 
(er =) So (Er 5) 
On en conclut par passage à la limite sur m : 
RE Se 
2(n+1)}?  n 7 2(n—1)? 


On peut alors écrire : 


1 1 1 1 1 
nd ibn I Ft < A ; 
Ps. ane 2\(n—1}? (n+1)? 
1 1 2 
n  2n? "| (n—1} 
On peut donc conclure que : 
1 1 
n  2n? 
est une bonne approximation de R,. 
5.5. Quelques sujets sur ce thème 
5.5.1. Exploration d’une suite (4). — On considère la suite de 
terme général : 
1 1 1 
h=gtat te 


Comparer expérimentalement S, et 7?/6. 
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Faire de même avec la suite de terme général : 


1 
Sn ir 
n 
5.5.2. Accélération de convergence. — On considère la suite de 
terme général : 
1 1 1 
Su — TE) + Les + — 


22 n2 


Comparer expérimentalement S, et 7?/6. 


Essayer de déterminer expérimentalement a de telle sorte que la suite 
de terme général : 


converge le plus vite possible. 


Le nombre a étant maintenant fixé à la valeur prévue précédemment, 
essayer de déterminer expérimentalement b pour que la suite de terme 
général : 


converge le plus vite possible. 


CHAPITRE 6 


DÉVELOPPEMENTS DÉCIMAUX DES 
NOMBRES RATIONNELS 


6.1. introduction 


Dans ce chapitre nous nous intéressons à des thèmes liés aux 
développements décimaux des nombres rationnels. Nous étudierons 
donc des fractions m/n qu'on supposera irréductibles et telles que 
0 < m < n. On trouvera une présentation détaillée du sujet dans l’article 
« Note sur les développements décimaux périodiques » (cf. [Rol)). 


6.2. La classification des comportements 


Nous rappelons que le comportement du développement décimal d’une 
fraction irréductible m/n est d’une des trois formes suivantes : 
— le développement décimal est fini, 
— le développement décimal est purement périodique, 
— le développement est mixte. 
Plus précisément : 
en est de la forme 2%5° auquel cas la fraction irréductible m/n 
admet un développement décimal fini : 


m 
— = 0.492: "5, 
n 


avec q; Z 0, ainsi qu’un développement décimal infini ne contenant 
que des 9 à partir d’un certain rang : 


m ! 
— = O.gige---g,999..., 
n 


où q! = qs — 1. 
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en n’est divisible ni par 2 ni par 5 auquel cas le développement 
décimal est illimité et purement périodique, c’est-à-dire que : 


m 
— 0.4192*-- 454102" Qs°°", 


ce qu’on note aussi : 

2 a 

— —= 0.q1q : --qs. 

n 
Nous adopterons la terminologie suivante : g142---q. est la partie 
périodique et s la période. 
e nest de la forme 2%50n, où n1 > 1 n’est divisible ni par 2 ni par 
5. Ce cas est un mélange des deux cas précédents. Il donne pour la 


fraction m/n un développement illimité périodique mixte, constitué 
d’une partie irrégulière et d’une partie périodique : 


fa > 
—. DUMAS; 


6.3. Les résultats classiques 


6.3.1. Cas d’un développement purement périodique. — Les 
décimales se calculent par des divisions euclidiennes successives ayant 
des restes non nuls : 

To — M, 


10ro = gin+r1 avec 0 <r1 < n, 


10r; = gn+r2 avec 0 < 79 < n, 
et par récurrence : 
10r; 1; = qn+r; avec 0 <7r; < n. 
Ce processus décrit parfaitement la division posée classique qu’on ap- 


prend à l’école primaire : ajout d’un zéro à la fin (multiplication par 10), 
recherche de g;, calcul du reste. 
Ceci fait jouer un rôle prépondérant aux deux suites qui interviennent 


dans les calculs précédents, c’est-à-dire d’une part la suite r des restes 
successifs : 


(SE (ri)i50 
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et d’autre part la suite des décimales : 
CE (dii>x ‘ 
Théorème 6.3.1. — Pour tout entier à > 0 on a : 


r; = 10%ro mod n. 


De cette proposition sur la suite r découle la proposition suivante sur 
les valeurs des termes de la suite q : 


Théorème 6.3.2. — Pour tout entier à > 1 on a : 


_ ut mod 2 
US n . 


Remarque: La proposition précédente nous permet de calculer une 
décimale particulière dont on précise le rang sans calculer les autres. Aïnsi 
on peut avoir une telle décimale même pour des grands nombres alors 
même qu'il n’est pas envisageable de les calculer toutes. Le comporte- 
ment du développement de la fraction irréductible m/n est indiqué par 
le théorème suivant : 


Théorème 6.3.3. — Les suites r et q sont périodiques de période s où 
s est le plus petit entier z > 1 tel que 10° = 1 (n). En particulier, la 
période s ne dépend pas de m (pourvu bien sûr que m soit premier avec 


La formule donnant la somme d’une progression géométrique nous per- 
met de reconstruire à partir d’un développement décimal périodique : 
0.g1Q2:--qs 
la fraction qui lui a donné naissance sous la forme : 


_ nds 
99...9 ? 


où le dénominateur comporte s chiffres 9. Bien entendu, à partir de là, 


O.g195:<"d3 


on se ramène à une fraction irréductible. 
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6.3.2. Cas d’un développement mixte. — Nous cherchons à 
connaître le comportement du développement décimal d’une fraction 
irréductible m/n où n est de la forme 25°n, avec n1 > 1 non divisible 
ni par 2 ni par 5 et a > 0 ou b > 0. 


Théorème 6.3.4. — Soit m/n une fraction irréductible telle que 0 < 
m<netn de la forme 2%5n; avec a > 0 ou b > 0 et n; non divisible 
par 2 ou par 5. Alors le développement décimal de m/n est périodique 
mixte c’est-à-dire de la forme : 


O.uiu2---UrG1Q2°*:qs, 


où ux £ qs. La partie uyu2 + - - ux est appelée partie irrégulière, sa longueur 
est k — max(a, b). La partie q1q2 - +: q, est la partie périodique, sa longueur 
s est le plus petit entier z > 1 tel que 10° =1 (rm). 


Comme dans le cas purement périodique, il est bien entendu possible 
de retrouver la fraction qui a donné naissance à un développement mixte : 


TE) 


———— 1 
Our uqiq "ds = 5x (arte ur + 99 ...9 


6.4. Approfondissements 


Dans toute cette section nous nous plaçons dans le cas d’un 
développement purement périodique. C'est-à-dire que nous étudions 
une fraction irréductible m/n avec 0 < m < n et n n'étant divisible ni 
par 2 ni par 5 (c’est-à-dire que n est un entier qui en écriture décimale 
se termine par 1,3,7 ou 9). 


6.4.1. La taille de la période. — Notons b(n) le nombre de nombres 
k premiers avec n tels que 1 < k < n. C’est aussi le nombre d'éléments in- 
versibles dans Z/nZ, c’est-à-dire encore le cardinal du groupe multiplica- 
tif des éléments inversibles de Z/nZ. On note (Z/nZ)* ce groupe. La fonc- 
tion @ est la fonction indicatrice d’Euler. On trouvera un développement 
complet de toute cette étude dans [Rol]. 


Théorème 6.4.1. — Soit a € Z/nZ*, alors aŸ() =1  (n). 
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Démonstration. — Si a est un élément de (Z/nZ)", il est facile de voir 
que l’application f de (Z/nZ)* dans lui même, qui à x associe f(x) = ax 
est une bijection. 

Soient 21,%2,:-: ,Té(n) tous les éléments de (Z/nZ)*. Comme f est une 
permutation on a successivement : 


g(n) d(n) 
Il: = II ads (0), 
i=1 i=1 


p(n) p(n) 
Il Tr; = a) El 2e CR), 
i=1 i=1 


et comme les éléments x; sont inversibles : 


a®) =1 (n) 
(m 
Théorème 6.4.2. — La taille s de la période est un diviseur de o(n). 
Démonstration. — Nous avons vu que cette taille s est indépendante du 


numérateur m, et plus précisément que c’est le plus petit entier z > 1 tel 
que 10*=1 (n). Faisons la division euclidienne de 6{n) par s : 


(n) = us + v, avec 0 < v < 8. 
On a alors : 
1090) = (10°)* x 10°, 
ce qui entraîne : 
10° =1 (n), 

et donc v = 0 (puisque v < 5). CI 

La taille s s'obtient en calculant toutes les décimales de la partie 
périodique, si cette taille est suffisamment petite. Mais lorsque n est 
grand on ne peut plus mener à terme le calcul en temps raisonnable. Si 


n est grand sans être trop grand (si n peut être factorisé) alors on peut 
calculer @(n), factoriser p(n) et tester les s possibles. C’est le cas si n à 
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une taille de 150 bits par exemple (de l’ordre de 50 chiffres décimaux). 
Cette démarche aboutit grâce à la décomposition de (n) : 


k 
o(n) = [[r°', 
=] 


et à l’algorithme suivant : 


s :— @{n); 
pour à := 1 jusqu’à k faire 
début 
tant que a; > let 10” =1 (n) faire 
début 
Si—S/p;; 
Q; := à; — l; 
fin ; 
fin ; 


retourner 5: 


Si n est vraiment trop grand et n’est pas d’une forme particulière, alors 
on ne peut pas le factoriser, on ne peut pas non plus calculer 6(n). Dans 
ce cas le problème du calcul de la période est hors de portée. C’est le cas 
actuellement si n à par exemple de l’ordre de 1024 bits ou plus. 


Dans le cas particulier où n est un nombre premier p alors d(p) = p—1, 
si bien que la longueur de la période est un diviseur de p — 1. Il arrive 
que cette longueur soit maximale, c’est-à-dire que s = p — 1. 


6.4.2. Fractions de même dénominateur. — On compare les 
développements décimaux purement périodiques de deux fractions 
irréductibles m1/n et m2/n où n n’est divisible ni par 2 ni par 5. On a 
vu que dans ce cas les deux périodes sont égales. Qu'en est-il des valeurs 
des parties périodiques ? Les résultats sont résumés dans le théorème 
suivant : 


Théorème 6.4.3. — S'il existe k tel que m2 = m1 X 10° (n), alors 
la partie périodique m2/n est à une permutation circulaire près la partie 
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périodique de m1/n. Plus précisément on obtient la partie périodique de 
m2/n en faisant tourner la partie périodique de m1/n de k positions sur 
la gauche. Sinon, les restes successifs de m2/n sont tous distincts de tous 
les restes successifs de mi/n. 


Quand la période est maximale, alors tous les restes possibles sont at- 
teints par m,/n et donc pour tout m2 premier avec n, la partie périodique 
de m2/n est à une permutation près celle de m1/n. C’est le cas par 
exemple de 1/7. Pour tout 1 < m2 < 7 la partie périodique de m2/7 est 
une permutation de celle de 1/7. Dans cette situation on peut dire que 
pour tout 0 < Mm1,Mm2 < n, il existe # tel que m2 = 10 X mu. 


1 
= — 0.142857 
“ 
2 
7 = 0.285714 2=1x102 (7) 
3 1 
7 0.428571 3=1%x 10 (7) 
4 
7 = 0.571428 ds 10 (0 
D 5 
70.714285 5=1%x 10 (gr) 
6 3 
70.857142 6=1%x 10 (7) 
6.4.3. Nombres de Sophie Germain. — Les nombres de Sophie 


Germains sont les nombres premiers q tels que p = 2 X q + 1 soit aussi 
premier. Pour le nombre p correspondant, il est alors très facile de calculer 
@(p) et de le factoriser, même si p est très grand. En effet : 


p(p)=p-1=2%X 4. 
En conséquence il est facile de voir que la période de 1/p ne peut être 
que : 
— 2 si et seulement si p = 11; 
— q si et seulement si 10%=1 (p); 
— p—1 si et seulement si 107% 1 (p) (et dans ce cas 101= —1 (p)). 
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6.4.4. Nombres de Mersenne généralisés. — On rappelle que les 

nombres de Mersenne sont les nombres de la forme 2° — 1. On peut 

considérer ces nombres comme cas particuliers des nombres de la forme : 
as —1 


a—] 


obtenus pour a = 2. Nous étudierons ici le cas des nombres obtenus pour 
a = 10, c’est-à-dire : 


10$ — 1 
5 
Un tel nombre s'écrit avec s chiffres 1 : 
10$ — 1 
= 11--.1. 
9 
Théorème 6.4.4. — Tout nombre n qui n’est divisible ni par 5 ni par 


2 divise un nombre dont tous les chiffres sont des 1. 


Démonstration. — Il suffit de développer 1/(9 X n) : 


1 —  Ujü Us 
= Quu2. Us = ——— , 
9xn 99 ...9 
ce qui prouve le résultat. ( 
6.4.5. La demi-période. — Soit p > 2 un nombre premier. On sup- 


pose que la période du développement décimal de la fraction 1/p est 
paire. Nous noterons 2s cette période. Ainsi : 


1 ——…—…—…—…—…———s 
— = 0.q1Q2: :-Q2s. 
P 


On constate alors que : 
s fois 
TR 
g192°*"Qs + GstiQe+2 * * * Gas = 99 ---9. 


Pour un nombre n qui n’est pas premier et tel que 1/n ait un 
développement purement périodique de période paire, ceci peut ne 
pas se produire, par exemple pour n — 21. Mais il arrive que ça se 
produise aussi : exemple n — 23 x 47 — 1081. En fait on a le théorème 
suivant : 
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Théorème 6.4.5. — Soit n un nombre entier qui n’est divisible ni par 
2 ni par 5 et dont la période du développement décimal est paire : 


1 ———_—— 
— = 0.g1g2:: "qe. 
n 


On a : 


19 * "se + Qsr1Qs42 * * Q2s — 99 -- 9 
si et seulement si : 
10° =-1 (n). 


6.5. Quelques calculs et essais 


Nous utilisons le logiciel xcas programmé en mode maple. 


6.5.1. Extraire les 2 et les 5 d’un dénominateur d. — Étant 
donnée une fraction n/d, dans la procédure qui suit on met cette frac- 
tion sous forme irréductible n’/q puis on « extrait les 2 et les 5 », du 
dénominateur q, c’est-à-dire qu’on calcule a et b tels que g = 2°x5°xm où 
m n’est divisible ni par 2 ni par 5. En outre on calcule aussi s — max(a, b), 
qui donne la longueur de la partie irrégulière du développement décimal 
de n/d. 


## Fonction val2s(n,d) 

## Entree : deux entiers positifs n et d (d>0) 

## Elle calcule q=-d/pgcd(n,d) au cas ou la fraction n/d ne serait 
## pas irreductible. 

## Soit q=2"a*5"b*m avec m non divisible par 2 ni par &5. 

#H# Sortie : les valeurs de a, b et s=max(a,b) dans une liste. 
val2s:=proc(n,d) 

local q,v,a,b,t; 

q:=d/gcd(n,d); 


t:={}; 
V:=d; 
a:=0; 


while irem(v,2)=0 
do 
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a:=a+i; 
v:=v/2 
od; 
t{[1]:=a; 
v:=d; 
b:=0; 
while irem(v,5)=0 
do 
b:=b+1; 
v:=v/5 
od; 
t[21 :=b; 
t[3]:=max(a,b) ; 
return t 
end; 


Voici des exemples d'exécution de cette fonction : 


val2s(1,6) 

1:04 
val2s(4,24) 

101 
val2s(1,140) 

22 
val2s(12,47) 

000 


6.5.2. Trouver la partie irrégulière d’un développement. — 

Étant donnée une fraction n/d, dans la procédure qui suit on calcule 

la partie irrégulière, si elle existe. C'est-à-dire que si le nombre s'écrit : 
n ss 2: 

à — Ouu2 - +: UxgiQ2""Qs; 


où u4 £ q,, alors la procédure renvoie : 


O.uiuo --- Uug. 
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Si la fraction n/d a un développement fini, c’est cette valeur qui est 


renvoyée. 


## Fonction irregu(n,d) 


## Entree : deux entiers positifs n et d avec n<d 
## Sortie : la partie irreguliere de n/d sous 


## forme decimale dans une liste 
irregu:=proc(n,d) 

local q,t,i,s; 

=}: 

q:=n/d; 

s:=val2s(n,d)[3]; 

for i from 1 to s 


do 
gq:=10%q; 
t[i] :=floor(q); 
q:=q-tli] 
od; 
return t 
end; 


Voici des exemples d'exécution de cette fonction : 


irregu (1,6) 


irregu(4,24) 


irregu(1,14) 


irregu(1,5) 


irregu(1,32) 


irregu(12,47) 
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Dans un deuxième temps, la fraction dont le développement est pure- 


ment périodique avec pour partie périodique celle de n/d est calculée. 


## Fonction fracper(n,d) 


## Entree : deux entiers positifs n et d avec n<d 


## Sortie : la partie périodique de n/d sous 


## forme de fraction dans une liste 


fracper:=proc(n,d) 


local t,a,b,s,x,u,p,q,l; 


1:={}; 
t:=val2s(n,d); 
a:=t[1]; 
bi=t{21;: 
s:=t[3]l; 
x:=10"s*xn/d; 
u:=x-floor(x); 
q:=d/(2"ax5"b) ; 
p:=u*q; 
1[1]:=p; 
1[21:=Q; 
return 1 
end; 

Exemples : 


fracper(12,47) 


H12 47] 
fracper(1,6) 
23 
fracper(1,140) 
57 
fracper(1,2) 
O1 


fracper(6313,49500) 
[53 99] 
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6.5.3. Fractions purement périodiques. — Ici on suppose que d 
n’est divisible ni par 2 ni par 5. On calcule la partie périodique du 
développement de n/d. Ceci ne s'applique que pour des d suffisamment 
petits. 


## Fonction period(n,d) 

## Entree : deux entiers positifs n et d avec n<d et 
## d non divisible ni par 2 ni par 5 

## sortie : la partie periodique dans une liste 


period:=proc(n,d) 
local t,a,r,i; 
t:={}; 

a:=n*x10; 
t[1]:=iquo(a,d); 
r:=irem(a,d); 


i:=2; 

while r<>n 
do 
a:=10%r; 


tlil:=iquo(a,d); 
r:=irem(a,d); 
i:=i+1 
od; 

return t 

end; 


Exemples : 


period(1,3) 
[3] 
period(1,7) 
142857 
period(1,11) 
(0 9] 
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period(2,11) 


H 8] 
period(1,13) 

076923 
period(4,13) 

307692 
period(5,13) 

384615 
period(6,13) 

461538 


6.6. Calcul d’une seule décimale 


Avec cette procédure on calcule une décimale isolée sans avoir besoin 
de calculer les autres. 
## Fonction decimale(k,n,d) 
## Entree : trois entiers positifs k, n et d 
## Sortie : la valeur de la k-ieme decimale 
## du developpement decimal de n/d 
decimale:=proc(k,n,d) 
local u,v,s,x,y,z,r,t,l; 
u:=n; 
v:=d\; 
if u>=v then 

u:=irem(u,v); 
fi; 
s:=val2s(u,v) [3]; 
if k<=s then 

x:=floor(10"kxu/v) ; 

y:=floor(x/10) ; 

z:=x-10*y 


t:=fracper(u, v); 
l:=period(t{[1],t[2]); 
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X:=k-s; 
y:=dim(1l); 
r:=irem(x,y); 
if r=0 then 
T:=y 
fi; 
zi=l[r] 
fi; 
return z 
end; 
Exemples : 
decimale(1,1,6) 
1 
decimale(3,6313,49500) 
7 
decimale(2007,12,47) 
0 
decimale(2007,12,1009) 
4 
decimale(2007,12,10019) 
6 
6.6.1. Nombres grands. — Voici des exemples de détermination de 


la période pour des nombres n relativement grands, mais pour lesquels 


on peut encore calculer et factoriser o(n). 


Nombre n : 


n = 256732889525456324556363824455622727 


Fonction d’Euler : 


@(n) = 171155158417286339647512680187384672 


Factorisation de phi(n) : 


2 x 3°? x 11° x 29 x 43? x 109 x 1427 x 11627153 x 50646958206521 
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Taille de la période : 
2 x 3x 11 x 29 x 43? x 109 x 1427 x 11627153 x 50646958206521 
— 324157497002436249332410379142774 


Nombre n : 
n = 4684296734077754585858421999 
99976543556735234875372852774373 
Fonction d’Euler : 
o(n) = 468429631828499565501852181 
874176326136092605087425689393456 


Factorisation de @(n) : 


2x 3x 7 x 13° x 17 x 37 x 89 x 149 x 1019 x 1289x 


805913 x 11050103 x 96096211 x 3813409399 x 230722618351 
Taille de la période : 
2x 3x7 x 13° x 17 x 37 x 89 x 149 x 1019 x 1289x 


805913 X 11050103 X 96096211 x 3813409399 x 230722618351 
— 585537039785624456877315227342 
72040767011575635928211174182 


= d(n)/8 


6.6.2. Nombres particuliers. — 


## Fonction neuf (p) 
## Entree : un nombre premier p>=7 
## Sortie : le nombre de 9 tel que pl9...9 
neuf :=proc(p) 
local 1; 
if p<7 then 
Disp("p<7 !") 
fi; 
l:=period(1,p); 
return dim(1l) 
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end; 
Exemples : 
neuf(13) 
6 
neuf(23) 
22 
neuf(127) 
42 
neuf(1277) 
638 
neuf(12739) 
4246 


6.7. Quelques sujets sur le thème des développements décimaux 
périodiques 


6.7.1. Développement décimal de période paire. — Le nombre 
1/7 a pour développement décimal 


1 
= = 0.142857. 
7 


Si on considère la première moitié de la partie périodique : 142, et qu’on 
l’ajoute à la deuxième moitié de la partie périodique 857, on obtient : 
999. Cela se produit-il dans d’autres cas (où la période est paire) ? 


6.7.2. Nombres de Mersenne généralisés. — Un nombre de Mer- 
senne généralisé de base a est un nombre de la forme 

Un 

a—1l” 


(a est un nombre > 1). Les nombres de Mersenne relèvent du cas où 
= 2, 


On va se placer dans le cas où a = 10. Un nombre de Mersenne de 
base 10 s’écrit sous la forme 11:--1 avec n chiffres 1. 
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Peut on trouver des nombres de Mersenne à base 10 qui soient des 
nombres premiers ? 


6.7.3. Nombres de Mersenne généralisés (bis). — Quels sont les 
nombres qui divisent un nombre de Mersenne généralisé de base 10 (un 
nombre de la forme 111---1)7? 


6.7.4. Nombres de Sophie Germain. — Peut-on trouver des 
nombres premiers q tels que p = 2q + 1 soit aussi un nombre premier 
(nombre de Sophie Germain) ? 

Dans ce cas quelle est la période du développement décimal de 1/p? 


6.7.5. Période maximale. — Existe-t-il des nombres premiers p tels 
que la période du développement décimal de 1/p soit égale à p—1? Dans 
un tel cas, soit 1 < k < p un nombre entier. Comparer le développement 
décimal de k/p et de 1/p. 
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